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Abstract. We present detailed modeling of the spectral energy distribution (SED) of the spiral galaxies NGC 
891, NGC 4013 and NGC 5907 in the far-infrared (FIR) and sub-millimeter (submm) wavelengths. The model 
takes into account the emission of the diffuse dust component, which is heated by the UV and optical radiation 
field produced by the stars, as well as the emission produced locally in star forming HII complexes. Radiative 
transfer simulations in the optical bands are used to constrain the stellar and dust geometrical parameters and 
the dust mass. We find that the submm emission predicted by our model cannot account for the observed fluxes 
at these wavelengths. Two scenarios that could account for the ‘missing’ submm flux are examined. In the first 
scenario dust additional to that derived from the optical wavelengths is embedded in the galaxy in the form of a 
thin disk. This additional dust disk, which is not detectable in the optical and which is associated with the young 
stellar population, gives rise to additional submm emission, and makes the total flux match the observed values. 
The other scenario examines the possibility that the average emissivity at submm wavelengths of the dust grains 
found both in a diffuse component and in denser environments (e.g. molecular gas clouds) is higher than the 
values widely used in Galactic environments. This enhanced emissivity reproduces the observed FIR and submm 
fluxes with the dust mass equal to that derived from the optical observations. In the second scenario, we treat the 
submm emissivity as a free parameter and calculate its nominal value by fitting our model to the observed SED. 
We find a dust emissivity which is ^ 3 times the often-used values for our Galaxy. Both scenarios can equally 
well reproduce the observed 850 pm surface brightness for all three galaxies. However, we argue that the scenario 
of having more dust embedded in a second disk is not supported by the near infrared observations. At 2.16 pm, 
the model images with a second dust disk reveal a prominent dust lane which is not present in the observations. 
Thus, the enhanced emissivity at submm wavelengths is a real possibility and the Galactic submillimeter dust 
emissivity may be underestimated. 

Key words. ISM: dust extinction — Galaxies: individual (NGG 891, NGC 4013, NGC 5907) — Galaxies: ISM — 
Galaxies: spiral — Infrared: galaxies — Submillimeter 


1. Introduction 

Edge-on spiral galaxies constitute a benchmark for multi¬ 
wavelength modeling of the dust content and the proper¬ 
ties of dust grains in spiral galaxies. Their unique orienta¬ 
tion allows for a direct detection of the dust distribution 
seen in absorption at optical wavelengths and in emis¬ 
sion in the far infrared (FIR) and submillimeter (submm) 
part of the spectrum. Realistic three dimensional radia¬ 
tive transfer (RT) modeling that was applied to a num¬ 
ber of nearby edge-on spiral galaxies (Xilouris et al. 1997; 
1998; 1999) has been able to determine the stellar and 
dust parameters that best fit the optical surface bright¬ 
ness of these galaxies. According to these simulations, a 
typical spiral galaxy contains ~ IO^Mq of dust grain ma¬ 
terial distributed in an exponential disk with a scaleheight 
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about half that of the stars and a scalelength about 1.5 
times the stellar one. Furthermore, the extinction law at 
optical and near infrared (NIR) wavelengths calculated 
for these galaxies matches very well the extinction law 
observed for our Galaxy, indicating common dust prop¬ 
erties among spiral galaxies. The optical thickness of the 
galactic disk, parameterized by the central face-on optical 
depth, has values of the order of unity in the B-band. 

Complementary to the optical modeling, studies of the 
FIR emission of edge-on spiral galaxies have also been car¬ 
ried out. In the studies of Popescu et al. (2000), Bianchi et 
al. (2000), Misiriotis et al. (2001), the FIR/submm energy 
output of a spiral galaxy is determined by equating it with 
the energy the dust grains absorb from the UV and optical 
photons. Popescu et al. (2000) found that their RT model 
submm fluxes were significantly lower than the observed 
values. In order to overcome this problem and explain the 
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FIR/submm spectral energy distribution (SED) of spirals, 
Popescu et al. (2000) used a second dust disk, with scale- 
height matching that of the young stellar population found 
in our Galaxy. However, recent studies question the valid¬ 
ity of the values that have been widely used so far for 
the FIR/submm emissivity of the dust grains. In partic¬ 
ular, a wide range of values (of an order of magnitude) 
are used for the submm emissivity (see Alton et al. 2004 
and Hughes et al. 1997 for reviews of the emissivity values 
found in the literature). In del Burgo et al. (2003), the 
FIR properties of dust in high-latitude regions are exam¬ 
ined using ISOPHOT maps and an enhancement (about 
four times) of the emissivity of the big grains with respect 
to those in the diffuse interstellar medium is indicated. 
In the study of Alton et al. (1998) and in a more recent 
study (Alton et al. 2004), the submm emissivity was cal¬ 
culated for the edge-on galaxies NGC 891, NGC 4013 and 
NGC 5907 by comparing the distribution of their visual 
optical depth with their 850 /rm emission. This analysis 
produced an emissivity at 850 /rm which is about four 
times the widely adopted value of Draine & Lee (1984). 
Alton et al. (2004) argue that since the submm emission 
closely follows the distribution of molecular gas, the rela¬ 
tively high emissivity values might be due to dust situated 
mainly in molecular gas clouds where the enhanced den¬ 
sity is conducive to formation of amorphous, fluffy grains. 
Such grains are expected to possess high emissivity val¬ 
ues. In a recent study of the large-scale variations of the 
dust optical properties in the Galaxy (Cambresy et al. 
2005), the authors suggest that this type of dust grain is 
more common than previously thought since they would 
be formed even at low extinction and not only in dense 
cold clouds. 

In this paper we use a three dimensional model of the 
emission of the dust grains in order to fit the FIR/submm 
SED of three edge-on spiral galaxies (NGG 891, NGG 
5907, and NGG 4013) already modeled in the optical wave¬ 
lengths by Xilouris et al. (1999). Two different scenarios 
concerning the dust properties are examined. In the first 
scenario (hereafter “1-disk” model) the dust is distributed 
in a single exponential disk (that derived by Xilouris et al. 
1999) with the FIR/submm emissivity treated as a free 
parameter so that the model SED matches the observed 
one. In the second case (hereafter “2-disk” model) an ad¬ 
ditional dust disk of smaller scaleheight with respect to 
the main dust disk is invoked. This disk is associated with 
the young stellar population and it contains, like the main 
disk, dust grain material with properties as described in 
Draine (2003). The sum of the two dust disks then gives 
rise to the FIR/submm flux that matches the observed 
SED. We show that both of these models are able to re¬ 
produce the observed surface brightness at 850 /rm. The 
“2-disk” model is tested for consistency by comparing the 
model K-band image with the observed one. The model 
K-band image exhibits a prominent dust lane which is not 
seen in the observations. From this test we conclude that 
the additional dust, if such is needed, cannot be in the 
form of a second dust disk. We cannot exclude the possi¬ 


bility that some dust is in the form of clumps. An upper 
limit for the amount of dust that can be located there has 
been placed by Misiriotis & Bianchi (2002). 

This paper is arranged as follows: In Sect. [2a descrip¬ 
tion of our model is given. In Sects. |2 and 0 the resulting 
SEDs of the “1-disk” and “2-disk” models are presented 
and evaluated. A comparison between the two models is 
made in Sect. |2 while possible sources of uncertainty are 
enumerated and discussed in Sect. 12 Finally, our work is 
summarized in Sect.|2 


2. Model 

2.1. General description 

The procedure that we follow to create a galactic 
FIR/submm spectrum is along the lines of Popescu et al. 
(2000) and is briefly presented here. Our first task is to 
find the radiation field in which the grains are immersed 
and from which they draw energy. For this reason we per¬ 
form RT calculations for every point in the galaxy and for 
all directions. These calculations need to be performed in 
all the wavelength regimes where light extinction by dust 
grains is taking place; namely, in the UV, optical and NIR. 
Gonsidering that the extinction properties of the dust in 
these regimes are well-known (Draine 2003), the RT cal¬ 
culations accurately provide us the power that the grains 
absorb per unit galactic volume. By adopting an appro¬ 
priate thermal emission law, equating the absorbed with 
the emitted power per unit volume, and assuming ther¬ 
mal equilibrium we are able to calculate the temperature 
distribution of the dust T(r, z), where r is the radial di¬ 
rection and z is the direction perpendicular to the galactic 
plane. Knowing the temperature distribution, we obtain 
the FIR/submm SED by integrating the emitted power 
Wem(A,T) over the entire galactic volume and converting 
it into an observable flux. The creation of such a model 
needs a priori assumptions of the way that the galactic 
stellar and dust components are distributed as well as the 
FIR/submm dust emission law. 


2.2. Stellar and dust distributions 


The geometry of the different components (stellar and 
dust) of the galaxies examined in this study is based on 
the analysis of Xilouris et al. (1999)| that simulates very 
accurately the optical images of these galaxies. The actual 
model that we use is briefly described below. 

The old stellar population is distributed in an expo¬ 


nential disk as well as a de Vaucouleurs (1953) bulge. In 
addition, a young stellar population is used as the main 
UV heating source of the dust grains. This component is 
distributed in a thin exponential disk in the galactic plane 
with the same scalelength as that of the old stellar popu¬ 
lation and a scaleheight of 90 pc as indicated by studies of 
the Galactic stellar population (Mihalas & Binney 1981). 
Quantitatively, the emission coefficients of a stellar disk 
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and a de Vaucouleurs (1953) bulge, rjs and rjb respectively, 
at a galactic position (r, z) and wavelength A are 


77s (A, r, 2 ) = 77s(A, 0,0) e e (1) 

r]b{X,r,z) = 776(A, 0,0)5“'^/® (2) 


In Eq. Vs and Zg give the scalelength and 

scaleheight of a stellar disk respectively. In Eq. m, 
B = + z'^{ajby(ve and depends on the effective ra¬ 

dius Ve and the major and minor semi-axes a, b of the 
bulge. 

The dust distribution within the galaxy is described 
by the extinction coefficient, which is thought to obey the 
exponential law 



fcext(A,r,0) = fcext(A,0,0) e-l"l/"T 


(3) 


The quantities and Zd are the scalelength and the 
scaleheight of the dust respectively. The central opti¬ 
cal depth of the model galaxy seen edge-on is given by 
r® = 2A:ext{A, 0,0)rd (Xilouris et al. 1997). In addition 
to this “main” dust disk, which is derived by modeling 
the optical appearance of edge-on galaxies (Xilouris et al. 
1999), in the case of our “2-disk” model, we use a sec¬ 
ond thinner dust disk of the same scalelength as that of 
the “main” dust disk and a scaleheight of 90 pc since this 
dust is assumed to be associated with the young stellar 
population (see above). 

The opacity (mass extinction coefficient), Kext(A) is de¬ 
fined as «;ext(A) = fcext(A,r, 0 )/pd(r, z), with pd{r,z) be¬ 
ing the density of the dust grain material inside the 
galaxy. It is connected to the mass absorption coefficient 
as Atabs(A) = (1 — a;)«;ext(A), where to is the scattering 
albedo. Taking into account that Kabs(A) = Kem(A) for 
A a where a is the radius of a typical dust grain (i.e. 
in the FIR and submm wavelengths), we derive that the 
mass emission coefficient (a grain’s emission cross section 
over its mass) is Kem(A) = (1 — a;)Kext(A). Additionally 
to the mass emission coefficient we can define the dust 
grain emissivity as Q(A) = (4/3)a/9«;em(A) (Whittet 1992), 
where p is the material density of the interstellar grains. 

We use the values of Xilouris et al. (1999) for 
all the quantities that can be inferred from an 
optical/NIR modeling of the galaxies. These are: 
77 s(A, 0, 0), 77 (,(A, 0,0), A:ext(A, 0,0) (for A in the optical/NIR 
waveband), Vg, Zs,re,a/b,rd, Zd, and the inclination angle 
of each galaxy. 

We parameterize the UV luminosity Ljjv of a spiral 
galaxy using the star formation rate (SFR) in 
For this, we use the formula of Kennicutt (1998) as de¬ 
scribed in Popescu et al. (2000)|| However, the UV lumi¬ 
nosity that heats the diffuse dust in a galaxy is less than 
that provided by the above equations for a given SFR. 
The reason is that a non-negligible fraction F of Luv can¬ 
not be radiated away from the young stars due to heavy 
obscuration of certain sightlines by the surrounding HII 
regions (Popescu et al. 2000). Thus, the fraction of Ljjv 
that is available to be absorbed by the diffuse dust is 1 —F. 


Fig. 1. SED modeling of NGC 4013 using different emis¬ 
sion laws. The solid and the dotted lines are derived with 
the emissivity described in Draine (2003) for Rv = 3.1 
and Rv = 5.5 respectively. The dashed-dotted line is 
the SED that we derive with the use of the model of 
Bianchi et al. (1999)[ and the short-dashed line is the SED 


produced with the Weingartner & Draine (2001) model 
for i?y=3.1. In all the cases we have set SFR = 1.8 
M 0 yr“^, and F = 0 , since we are interested only in the 
behavior of the diffuse dust. 



Fig. 2. The Draine (2003) emission coefficient times A^ 
for Rv = 3.1. The curve is nearly flat in the regime of our 
concern (50 < A < 1500/rm). 


We add the HII region contribution to the SED of the dif¬ 
fuse dust by using the spectrum of the Galactic region 
G45.12-1-0.13 (Ghini et al. 1986) as a template and nor¬ 
malizing its luminosity to FLjjv- 

2.3. FIR/submm emission law 

Having identified all the ingredients of our model and for¬ 
mulated their distributions, we still have to assume a dust 
FIR/submm emission law in order to find the SED, as 
mentioned in Sect. 2.1. 
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X (|rm) 

Fig. 3. Same as Fig. 1. The dashed-dotted line is the SED 
derived from our model with SFR = 1.8 M0yr“^, F = 0, 
and Co = 987. The solid line is identical to that of Fig. 1. 


The selection of the appropriate values for the dust 
emission coefficient is a major decision in our modeling. 
As we will soon demonstrate, all the emission coefficients 
that have been found so far for the diffuse dust of the 
Milky Way fail to explain the FIR/submm SED of the 
galaxies in our sample, if we assume that the amount of 
dust derived from the optical images is correct. This leads 
us to investigate other ways that can produce the observed 
SED. 

We begin by using the values that Draine (2003)| pro- 
poses for the Milky Way. These values were calculated 
by taking into account a mixture of graphites, silicates, 
and polycyclic aromatic hydrocarbon molecules (PAHs), 
as well as a distribution of grain sizes (Weingartner & 
Draine 2001). In particular, we use the tabulated, more de¬ 
tailed data of the website cited in Draine (2003) Several 
values of Kem found for different grain abundances rela¬ 
tive to H are available and described as a function of Ry- 
According to Draine, the quantity Ry = AyI{A b — Ay) 
shows how steep the slope of the extinction Ax is in the 
optical waveband. The extinction is expressed in terms 
of the flux Fx as = 2.5 log^jj(F°l>served/^u„exti„guished^^ 

We use the cases Ry = 3.1 and Ry = 5.5. The Ry = 3.1 
model is considered representative of the mean Galactic 
obscuration (Whittet 1992), while the Ry = 5.5 case 
roughly agrees with the mid-infrared extinction law at ~5 
/rm, as observed by Lutz et al. (1996)| in the Galactic cen¬ 
ter direction. 

The emitted power per galactic unit volume and A, 
Wem, is equal to the flux of the (gray-body) radiation 
through each grain’s spherical surface, times the emitting 
surface, times the number density of the grains. When we 
combine the above with the definition of the mass emission 
coefficient we obtain the following emission law 


Wem(A,r) = inpd Kem(A) B{X,T), (4) 

where B{X,T) is the Planck function and T is the tem¬ 
perature of the dust. As already mentioned, we equate the 


integral of Eq. (4) to the absorbed energy per unit volume 
to find the temperature of the dust grains. Then, having 
the T distribution in the galaxy, we return to Eq. (4) to 
calculate the emitted power as a function of A. 

In Fig. n we show as a solid line the SED for NGG 
4013 produced with Draine’s (2003) emissivity, SFR = 1.8 
M0yr“^, F = 0, and Ry = 3.1. The value of SFR was 
chosen so that we obtain a good fit to the 60 and 100 p,m 
data (see Sect. 3). The dotted line is the same as the solid 
one but for SFR = 1.8 M0yr“^, F = 0, and Ry = 5.5. 
The reason we set F = 0 at this point is to compare 
several emission models that apply only to the diffuse dust. 
Any conclusion for the submm emissivity derived from this 
comparison will not be changed with the addition of the 
HII regions’ emission, because it mainly contributes at the 
FIR wavelengths (see Popescu et al. 2000). Glearly, if we 
have included in our modeling all the dust that exists in 
NGG 4013, these emissivity models cannot account for the 
galaxy’s submm SED. 

We reach the same conclusion when we use emissiv¬ 
ity laws found by other authors. The dashed-dotted line 
in Fig. n is the SED that we derive with the use of 
the model of Bianchi et al. (1999) This model uses larger 
(and therefore colder) grains than Draine (2003), SFR = 
1.8 and F = 0. The dashed line is the SED pro¬ 

duced with the Weingartner & Draine (2001) model for 
i?y=3.1, SFR = 1.8 M0yr“^, and F = 0. In all cases, 
the observed flux at 450 and 850 ^m is significantly higher 
than the above models. Similar results were reported by 
Popescu et al. (2000)||and Misiriotis et al. (2001)| 


The approach of Bianchi et al. (1999)| in d’^ving the 
emissivity is conceptually simpler than that of the other 
authors, but still trustworthy. Bianchi et al. (1999) used 
a model that is independent of grain size and composi¬ 
tion and that does not take into account the PAH emis¬ 
sion in order to explain the FIR flux maps of the Galactic 
hemispheres. The wavelength dependence of the emission 
coefficient, as inferred from that study, has the form 


Kext(f^) M00^mN2 
RAJ - 760 V A 7 ’ 


(5) 


with /vext(F) being the V band opacity. 

Driven by the fact that even the most sophisti¬ 
cated models do not differ much from a 1/A^ law in 
the range 50 < A < 1500/rm (see Fig. |21l, we adopt the 


Bianchi et al. (1999) assumptions of grains of an average 


size and composition, that are gray-body emitters. We 
consider the emission coefficient to be 


Kext(lA) ^ 100^m\/3 
Gn I A ^ ’ 


( 6 ) 


Co V A 

where /3 = 2 , Gq is a parameter to be determined, and 

^ (Draine 2003, Ry = 3.1). 
1-disk” model are SFR, F, 


cm^ gr 


Kext(E) = 2.52 X 10^ 

Thus, the parameters in our 
and Co- 

These assumptions lead to the gray-body emission 
form that we adopt, namely 


(A,r) = b{X,T), 


V A 


(7) 
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Table 1. The best fitting parameters for our “1-disk” 
model. 


Table 3. The best fitting parameters of our “2-disk” 
model. 


galaxy 

SFR (Moyr-^) 

F 

Co 

NGC 891 

5.80 

0.28 

330 

NGG 4013 

1.80 

0.17 

187 

NGG 5907 

4.10 

0.16 

236 


Table 2. Ratios xb,xw, 

and 

Xd of our emis- 

sion coefficient to that 

of 

Bianchi et al. (1999)| 

Weingartner & Draine (2001) 

and 

Draine (2003)1 re- 


galaxy 

SFR 

(Moyr-i) 

F 

Mdi 

(Me) 

Md2 

{Me) 

NGG 891 

3.40 

0.40 

5.6x10' 

1.3x10“ 

NGC 4013 

0.80 

0.32 

4.5x10® 

1.6x10’'' 

NGC 5907 

2.40 

0.26 

1.5x10''' 

3.5x10'^ 


3. Results of the “1-disk” model 


spectively, at 850 pm. 


galaxy 

Xb 

Xw 

Xd 

NGC 891 

2.3 

2.1 

2.8 

NGC 4013 

4.1 

3.7 

4.9 

NGC 5907 

3.4 

3.0 

3.9 


with /3 = 2. In Sect. 6 we will examine other values of [3 
to see how sensitive our conclusions are to such changes. 
Furthermore, we will comment there on the effects of the 
PAHs and the transiently heated, very small grains (VSG), 
which we did not take into account in our modeling. 

What SED would our model give for NGC 4013 in 
comparison to those shown in Fig. 1? In Fig. 3 we show as 
a dashed-dotted line the SED produced with our model 
for SFR = 1.8 M0yr-\ F = 0, and Co = 987. It is 
remarkably close to the solid line, which is identical to 
the solid line of Eig. 1. This means that a simple model 
with the appropriate Kem (i-e. appropriate grain size and 
composition) can mimic a more complicated model and 
justifies our selection. 


2.4. The ‘‘2-disk” model 


Eor our “2-disk” model, we keep the stellar disks (old and 
young) the same as in our “1-disk” model. The dust, how¬ 
ever, is now distributed in two disks. One, with mass M^i, 
is identical to that of our “1-disk” model. The second, with 
mass Mii 2 to be determined, has scalelength equal to 
that of the first disk. Its scaleheight Zd is taken equal to 90 
pc, as in Popescu et al. (2000) As these authors explain, 


the second dust disk must have a small scaleheight to be 
undetected in the optical bands. This is true, at least for 
wavelengths less than ^ 2pm, because the “main” dust 
disk highly obscures any other feature (the second dust 
disk included) which is close to the galactic plane. The 
specific value of 90 pc was used by Popescu et al. (2000) 
because they associated the second dust disk with the 
molecular cloud distribution in the Milky way and pro¬ 
posed that it has to be connected with the young stel¬ 
lar population. The addition of the second disk obviously 
leads to an increase of the central extinction coefficient of 
Eq. 0. Another difference to the “1-disk” model is that 
the FIR/submm emissivity is taken from Draine (2003). 
The free parameters of the “2-disk” model are then SFR, 
F, and Md 2 - 


Following the procedure described in Sect. [3 we deter¬ 
mined the parameters that best fit the observational data 
in the FIR/submm waveband. Their values are given in 
Table 1 for the three galaxies of our sample. 

In Fig. 0 we show our computed SEDs that best fit 
the observational data. Eor NGC 891 the observed fluxes 
are taken from Pupae et al. (2003a)[ with the exception 
of the ISO and the IRAM measurements which are taken 
from Popescu et al. (2004) and Guelin et al. (199^ re¬ 
spectively. Eor NGC 4013 and NGC 5907 the IRAS data 
are taken from Soifer et al. (1989) while the SCUBA data 
are presented in Alton et al. (2004). The 1.2 mm flux of 
NGC 5907 is reported in Dumke et al. (1997). The short- 
dashed lines give the contribution to the SEDs of the HII 
regions, the dashed-dotted lines give the contribution of 
the diffuse dust and the solid lines represent the total com¬ 
puted SEDs. It is evident that our model SEDs fit very 
well the observational data. 

It is also evident from Table 1 that the average 
value of our parameter Cq is 251 and corresponds to 
Kem(850/im) = 1.4 cm^ when using Kext{V) = 

2.52 X 10"^ cm^ gr”^, as the Draine 2003, Ry = 3.1 model 
prescribes. In Table 2, we compare for each galaxy the 
850 /im emission coefficient that we find (from Eq. |3) to 
values used in Galactic emission studies. More specifically, 
we present the ratios x b , Xw, and xd oi our 850 fim value 
to those of Bianchi et al. (1999), Weingartner & Draine 
(2001), and Draine (2003) respectively. The average of the 
above ratios is 3.4, but there is a spread of values due to 
the assumptions of the Galactic models and to the vari¬ 
ability of the measured extinction amongst the Galactic 
sightlines. 

Where does our result for Kem(850/rm) stand relative to 
those for other environments? The average of the values of 
ccb in Table 13 is comparable to Bianchi et al. (2003) when 
they studied Barnard 68, a dark cloud in the foreground 
of the Galactic bulge. They concluded that the most 
characteristic «:em(850^m)/Kem(U) ratio is 4.0 x 10“^ in¬ 
stead of 1.8 X 10“® that results from Eq. 0. The study 
by [del Burgo et al. (2003) of high-latitude Galactic inter¬ 
stellar regions at 200 /im led to similar findings. These 
authors suggested that the large and cold grains (simi¬ 
lar in size to the ones we use and of temperature ~ 14 
K) have an increased FIR emissivity by a factor > 4. 


Ossenkopf & Henning (1994) studied the dust emissivity 


in protostellar cores where grains are often covered with 
a mantle of ice. They found that at 850 /im the emis- 
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Fig. 4. The observational data of the SED of NGC 891, NGC 4013, and NGC 5907 and our “1-disk” model fits (solid 
lines) to them. The short-dashed and dashed-dotted lines give the contributions of the HII regions and the diffuse dust 
respectively. 



Fig. 5. Same as in Fig.Q] but for our “2-disk” model. 


sivity is about 4-5 times higher than that of the in¬ 
terstellar medium as given by Draine & Lee (1984). In 
addition to this they presented a less realistic scenario 
of cores with bare grains (without ice coating), where 
Kern = 3.5 cm^ gr ^ at 850 /rm, for atomic hydrogen den¬ 
sity TiH = 10® cm“®. This is an order of magnitude higher 
than Draine’s (2003) value (/tem = 0.382 cm^ gr ^). The 
crucial question, though, is to what degree emissivities 
derived from such environments can be used for a spiral 
galaxy as a whole. 


Other authors (Dumke et al. 1997) believe that X has 
so many uncertainties on its own, that it is safer to cal¬ 
ibrate its value with the aid of the dust emission. More 
details on the accuracy of the X factor can be found 
in Maloney &: Black (1988)|and|Arimoto &: Sofue (1996) 


Nevertheless, [James et al. (2002)| claim that the uncer- 
tainty attached to their result is a factor of 2, in which 


fjames et al. (2002) calibrated the emission coefficient 
at 850 /im for a sample that included several types of 
galaxies and found that Kem = 0.7 cm^ gr“^. The method 
applied was to express the dust mass in two ways, one 
in terms of the gas mass and the metallicity, the other 
in terms of the gray-body emission law, and then to 
equate the two expressions. Thus, they managed to solve 
for the emission coefficient independently of the dust 
mass. Although this is the theoretically optimal treat¬ 
ment, their method assumes that the amount of met¬ 
als locked up in dust grains is constant for all galaxies. 
They also rely on the CO to H2 conversion factor X. 


case it is closer to our result. Dunne et al. (2000) cal¬ 
culated the dust masses for the Scuba Local Universe 
Galaxy Survey (SLUGS) as suming Kem = 0-77 c m^gr~^, a 
value very similar to that of James et al. (2002) Recently, 
Seaquist (2004) found that for a revised conversion factor 
A, the masses of the SLUGS galaxies have to be reduced 
by 25-38%, which would lead to a similar increase in their 
value of Kem, reducing the difference from our result. 

Our calculations for the values of Kem (850 Atni) can be 
directly compared with the studies of jAlton et al. (1998) 
and Alton et al. (2004) where the three galaxies studied 
here are also included. Although their method is based 
on simple assumptions (they compare the model derived 
optical depth with the 850 /im flux density), they derive 
a value for Kem(850/im) which is 4 times that of Draine 
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Fig. 6. The “1-disk” and “2-disk” theoretical images (of 
NGC 5907) at 850 /rm and their percentage difference are 
displayed from top to bottom. The deviation is illustrated 
by the key at the bottom. It is obvious that in the submm 
the two model images are very similar to each other. 


& Lee (1984). This result comes in good agreement with 
the conclusions of our method, which uses a much more 
realistic way of calculating the submm emission of the 
galaxy (by computing the temperature distribution in a 
self consistent way for every point inside the galaxy, in 
contrast with the method of Alton et al. 1998 and Alton 
et al. 2004, where they use simple grey-body fits to derive 
the dust temperature). 

Another recent estimate of the 850 ^m emission co¬ 
efficient in nearby spirals is given in the recent study of 
on the face-on galaxy M51. These 
/rm map of M 51 and estimate that 
«;em(850 ^m) is 1.2 cm^ gr“^ (about two times higher com¬ 
pared with values widely used in the literature) based on 
a canonical gas-to-dust ratio. 

4. Results of the “2-disk” model 

Once more, we followed the prescriptions of Sect. [3 to 
determine the parameters that give the optimal fit to the 
observational data. Their values for this “2-disk” case are 
presented in Table 3 and their computed SEDs are shown 
in Fig.|3as solid lines. As in the “1-disk” model, the SEDs 
fit the observational data very well and do not distinguish 
between the two possibilities. 

Erom Table 3 we see that we need a total dust mass 
Md which is 3.6 times that of the “1-disk” model (on aver¬ 
age for the three galaxies), while in Table 2, an enhanced 
emissivity of about 3.9 times (again on average) that of 
Draine’s model is needed in order to account for the ob¬ 
served flux. This is not an unexpected result. Under the 
simplifying assumption that the FIR/submm SED can be 


Meijerink et al. (2005) 
authors present an 850 


described by a gray-body law, the quantities Md and Kem 
have an inverse proportionality relation for a given 850 /im 
flux (e.g. Dunne & Eales 2001). In reality, this assumption 
is only valid for optically thin environments. The optically 
thicker the environment gets, the more the results deviate 
from what the above reasoning indicates. 


Popescu et al. (2004) found that the second dust disk 
of NGC 891 contains mass equal to 7 x 10^ Mq and thus a 
total dust mass 2 times the value of Xilouris et al. (1999) 
(instead of 3.4 that we found). This difference may 
originate from the adopted emission law. Popescu et 
al. (2004) described the grain emission following the 


Laor & Draine (1993) recipe while, as already mentioned, 
we selected the more recent values of Draine (2003) 


Except for the total dust mass, a notable parameter 
now is the SFR, which is approximately half of that needed 
in the “1-disk” case to provide the required UV radiation 
to heat the dust. This is reasonable because as the num¬ 
ber of emitting particles increases, so does their emitted 
power, if they are in an optically thin or moderate envi¬ 
ronment. 


5. Comparison between the two models. 

By either adjusting the dust emissivity or the dust mass 
we succeed in fitting the FIR/submm SED. In this Section 
we compare the two models. 

We start with the submm regime. In Fig. |3we show 
the 850 /xm theoretical images of the “1-disk” model (top), 
the “2-disk” model (middle) and the percent difference of 
the second to the first (bottom). The fractional deviation 
between the two images is less than 20% throughout the 
galaxy. Unfortunately, the data is inconclusive because of 
the poor SCUBA resolution and because of the short spa¬ 
tial extent of the galaxies along the vertical direction. 

Nevertheless, the data resolution is sufficient for the 
comparison of the 850 /xm flux radial distribution, along 
the major axis, with that predicted by our models. Since 
we know the three dimensional distribution of the emit¬ 
ted power within the galaxy, the radial profiles are easily 
created by projecting all galactic positions’ flux to the 
central plane. Then, the models are smoothed to the data 
resolution. The results are given in Fig. [7| (left panels). 
The solid lines correspond to the “1-disk” model while the 
dashed ones to the “2-disk” model. Both models provide 
acceptable fits if one neglects the peculiarities of each in¬ 
dividual galaxy which cannot be reproduced by a smooth 
model. In the same Figure, the temperature distribution 
of the grains is given. Both cases agree with other authors 
(Dumke et al. 1997; Dunne et al. 2000; Dupac et al. 2003a; 
Dupac et al. 2003b; Alton et al. 2004). 

We compare the NIR appearance of one of the galaxies 
in our sample (namely NGC 891) when the two different 
models (“1-” and “2-disk”) are considered. The extinction 
caused by the dust in these wavelengths, and especially in 
the K-band, is low, making any additional dust features 
easily detected. For NGC 891, the central edge-on optical 
depth in the K-band is t](.=2.7 (as derived from Xilouris 
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radial distance (kpc) t (kpc) 


Fig. 7. Flux radial profiles and temperature radial distribution for NGC 891, NGC 4013, and NGC5907. The solid and 
the dashed lines correspond to the “1-disk” and the “2-disk” models respectively. The temperature profiles of both 


models are plotted for several galactic heights, namely, z =i 
pc for NGG 891 and NGC 5907 from top to bottom. 

et al. 1999). Thus, a very weak dust lane is seen in the 
K-band image (see the upper panel of our Fig. EJ. The 
“1-disk” model fits very well the K-band image of NGC 
891 (Xilouris et al. 1998; see also the middle panel of our 


', 400, and 1600 pc for NGC 4013 and 0 =0, 600, and 2400 


Fig.lHJ. The “2-disk” model on the other hand has a cen¬ 
tral, edge-on optical depth equal to 15.6 in the K-band and 
creates a very prominent dust lane (Fig. El lower panel). 
Such a high value for tk comes from the fact that we have 
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added the contribution of the second, thin, dust disk to 
that of the first disk to find the total optical depth. Since 
the second disk has a scaleheight of 90 pc, it is ^ 3 times 
thinner than the first one (Xilouris et al. 1999), contribut¬ 
ing significantly to rf,. 

A more direct comparison between the two different 
models (“1-disk” and “2-disk”) is presented in Fig. 0 In 
this plot we show the vertical profile of the K-band image 
of the galaxy averaged over a region of 430 arcsec along 
its major axis (stars) together with the vertical profiles 
of the “1-disk” model (solid line) and the “2-disk” model 
(dashed line). The region where the profiles are averaged 
covers most of the galaxy detected in the K-band. This 
comparison shows the excellent fit of the “1-disk” model 
(solid line) to the observations; when a second dust disk is 
included, the dust content is overestimated giving a ver¬ 
tical profile that does not match the data. It is therefore 
impossible to hide a second dust disk which has compa¬ 
rable or larger dust mass than the first, even if it does 
not appear in the optical wavebands. Of course, smaller 
amounts of dust may have gone undetected by our mod¬ 
eling. 

The K-band image modeling clearly favors the “1-disk” 
model. Although we only have K-band data for NGC 891, 
we have no reason to believe that this is a special case. 
Thus, our work indicates that the dust emissivity in spiral 
galaxies at FIR/submm wavelengths seems to be about 3 
times what is thought appropriate for the Milky Way. 


6. Possible sources of uncertainty 


It is important to check whether inaccurate knowledge of 
certain quantities can lead to different conclusions. 

So far we have only considered the possibility that any 
additional amount of dust (to that in the first disk) lies in 
a diffuse state (i.e. in a second disk). We also have to inves¬ 
tigate the possibility of keeping only the first disk and dis¬ 
tributing some extra dust in dense and quiescent clumps. 
This scenario could be realistic because the clumps do not 
appear in optical images when not in a high number den¬ 
sity. In addition, due to their high concentration of grains, 
they may contain large amounts of dust and, thus, they 
may be able to account for the unobserved flux. We will 
investigate this scenario for clumps of different masses (of 
order 10^ — 10'^ Mq and of order > 10® Mq respectively). 

We argue that it is not reasonable to distribute dust in 
clouds similar to the dark cloud D of M17 (as denoted by 
Dupac et al. 2002), which belongs to the lower mass cat¬ 
egory. The reason is that the number of clouds needed to 
reproduce the observed SED renders the galaxy optically 
thick. In the case of NGC 4013, the emission of the second 
dust disk is 2.5x 10^^ erg s“^. Presuming that this flux now 
originates from quiescent clumps, we calculate their num¬ 
ber as follows. We fit a gray body of /3 = 1.9 and T=14 K 
to the data of Dupac et al. (200^ as these authors pre¬ 
scribe. This allows us to find the FIR/submm emitted 
power of the M17 cloud D, which is equal to 2.8x10®® 
erg s” 



Fig. 8. K-band image of NGC 891 (upper panel), “1-disk” 
model (middle panel), and “2-disk” model (lower panel). 
The predicted images have been smoothed according to 
the seeing conditions of the dataset. 



Fig. 9. The vertical profile of the K-band image of the 
galaxy averaged over a region of 430 arcsec along its major 
axis (stars). This region covers most of the galaxy detected 
in the K-band. Along with the data we present vertical 
profiles of the “1-disk” model (solid line) and the “2-disk” 
model (dashed line). 
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Fig. 10. The left panel is the SED of NGC 4013 for the case of a /3=1.9 emission law (dashed-dotted line) and for the 
case of /3=1.5 (short-dashed line). For the sake of clarity, we plot in this panel the total emission (diffuse dust and HII 
regions emission together) for each case. The observational data are the same as those in Fig. 0] The right panel is 
the SED of the same galaxy plotted together with data corrected for VSG and CO emission. This panel is similar to 
those in Fig. 4 (HII regions, diffuse dust and total emission are plotted in short-dashed, dashed-dotted and solid line 
respectively). 


clouds to produce the observed SED. In order to check the 
effect of such a number of clouds on the optical thickness 
of NGC 4013, we need a rough estimate of the radius of 
the M17 cloud D. From Pupae et al. (200^ it can be in¬ 
ferred to be approximately 10" or 6.4 pc. When 8.9x10® 
clouds of this radius are uniformly spread on the galac¬ 
tic plane (one next to the other), they form an optically 
thick disk of radius 6 kpc. For NGC 4013 this corresponds 
to ~ 3rs, extending almost to the edge of the (baryonic 
component of the) galaxy and, thus, renders the galaxy 
opaque. Similarly to the diffuse case, the addition of dust 
has to be so high that it will inevitably appear in the 
optical or NIR images. 


Clumps of larger size and mass (probably associated 
with Giant Molecular Clouds) and their effects on the op¬ 
tical thickness of spiral galaxies have been treated in a 
more sophisticated manner by Misiriotis & Bianchi (2002). 
These authors showed that clumping can lead to an un¬ 
derestimate of the dust mass by 40% at most. In Sect. 0] 
we showed that the mass addition necessary to reproduce 
the observed SED is ~ 300%. Thus, these clumps only 
contain a very small fraction of the required dust mass. 
On the other hand, due to their lower temperatures, they 
will primarily emit in the submm. Thus, the question is 
to what degree will their contribution lower the value of 
Kem that we find for the diffuse dust. To answer that, we 
intend to implement radiative transfer models that take 
into account dumpiness. 


A second source of uncertainty is the value of the quan¬ 
tity (3 in the emission law (Eq. j^]). The average temper¬ 
ature of the dust in the galaxies of our sample is 16.5 K 
in the “l-disk” model. According to |Dupac et al. (2003b)| 
the appropriate value of /3 for this temperature is 1.9 
rather than 2. If this is the case, the optimal Kem of NGC 
4013 corresponds to Cq = 236, while the value of the pa¬ 
rameters SFR — 1.77 MQyr“^ and F = 0.17 remain the 
same. This results in ratios xb = 4.0 and xd = 4.8, which 


differ very little from the values given in Table [3 In this 
case, the total SED for NGC 4013 (the SED of both the 
diffuse dust and the HII region components) is given in 
Fig. ^1 (left panel, dashed-dotted line). 

Now we turn our attention to the emissivity at 100 ^m. 
Our average best Ht parameters for the (3 = 2 law result 
in a very high value of the emission coefficient at 100 ^m. 
Indeed, in that regime the dust emission is considered to 
be more accurately described by a /3 = 1.5 law. For this 
reason we also use Eq. Q with this new value for (3 and 
we find that the best fit is depicted in the values Gq = 645, 
SFR = 1.70 M0yr“^ and F = 0.13. The ratios xb and 
Xd are now equal to 3.4 and 4.1 respectively. The best fit 
is again given in the left panel of Fig. irni (short -dashed 
line). Our conclusion is that changing the value of j3 to 
1.5 does not have a significant impact on our results. 

Since both values of (3 for our FIR and submm bound¬ 
ary conditions (100 /rm and 850 /rm respectively) lead to 
a similar value of Kem(850 i^m), our results are proven to 
be robust. We are also conhdent that any other emission 
law of variable (3 that complies with these boundary con¬ 
ditions (like that of Reach et al. 1995) will lead to similar 
findings for the 850 /rm emission coefficient. We still prefer 
to use a constant (3 model because it has the advantage of 
being simple. 


Two other concerns are the possible contamination of 
the FIR data from molecular line emission and the effects 
of VSG. We start with the J = 3 — 2 emission since 
it contributes to the 850 /rm flux, from which we derive 
our conclusions. ||Papadopoulos &: Allen (2000) found that 
it accounts for 40% of the total emission in the starburst 
galaxy NGG 7469. Of course, it has to be much smaller for 
a quiescent spiral. Dumke et al. (2001) mapped the ^^GO 
J = 3 — 2 emission for different types of nearby galax¬ 
ies and found that the total emitted power in this tran¬ 
sition is 2.5 X 10®® erg s“® for NGC 891. This accounts 
for 1.2% of the flux we used at 850 /rm, but it is only a 
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lower limit because the ^^COJ = 3 — 2 mapping did not 
cover the whole area of emission. Dumke et al. (2001) ar¬ 


gue that the unobserved flux is less than a few percent of 
the total. Meijerink et al. (2005) who studied the submm 
emission in M51 discussed that the molecular contamina¬ 
tion at 850 /rm is at most 35% in the spiral arms and ap¬ 
proximately (less than) 12% for a homogeneous dust disk 
in the inter-arm regions. We use the average value, 24%, 
which is close to that of Alton et al. (2004) (20%). At 450 
/rm, Papadopoulos & Allen (2000) claim that the molec¬ 
ular contribution that comes from the J = 6 — 5 

line represents less than 2% of the total flux in a starburst 
environment. On the other hand, Alton et al. (2004) give 


an approximate correction of 10% at 450 /rm, which is the 
one we use. The correction of the ISO data for the VSG 
emission comes from the same source (Alton et al. 2004), 
is 62% at 60 /rm and 14% at 100 /rm, and is believed to 
be somewhat overestimated. Conclusively, we correct the 
60, 100, 450, and 850 /rm fluxes for our galaxies by reduc¬ 
ing them by 62%, 14%, 10%, and 24% respectively. The 
best fit for NGC 4013 is given in Fig. II01 right panel. The 
parameters now are SFR — 1.6 MQyr“^, F = 0.02, and 
Co = 262. The value of Cq is robust against contamination 
of our data, and again, our main results are unaffected. 


7. Summary and conclusions 

We presented detailed modeling of the FIR/submm SED 
of three edge-on spiral galaxies (NGC 891, NGC 4013, and 
NGC 5907). We demonstrated that the observed SED can 
be fitted equally well by a “l-disk” model, where the dust 
(having a submm emissivity ~ 3 times the value widely 
used in the Galaxy) is distributed in a single exponential 
disk, or by a “2-disk” model where additional dust (~ 2 
times more) is distributed in a second thinner disk associ¬ 
ated with the young stellar population in spiral galaxies. 

The FIR/submm spectrum of spiral galaxies was not 
able to distinguish between the two models. We showed 
however, that the “2-disk” model is unrealistic due to its 
intense dust lane in the K band, which is not present in 
the observations. 

Future research on the dust distribution within galax¬ 
ies (e.g. including more realistic cases like spiral arms, 
dust clumps, etc.) has to be conducted on both edge-on 
and face-on systems in order to investigate the existence 
of more dust hidden within the galaxies and, thus, obtain 
a better estimate of the properties of the dust grains. 
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